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I. INTRODUCTION 


The description of high-energy scattering processes depends on what is measured in ex¬ 
periments. In the inclusive scattering in which physical observables do not depend on the 
identity of the final-state particles, we sum over all the possible configurations of the particles 
at the parton level, or at the hadron level. If we consider the exclusive scattering in which 
the observables depend on the final-state hadrons, the hadronization process through the 
fragmentation functions should be implemented. In all the processes, factorization theorems 
constitute an important ingredient to separate nonperturbative effects from hard scattering. 

The factorization proof of high-energy scattering processes has been of great theoretical 
interest, and the factorization in full QCD was considered in the pioneering work of Collins, 
Soper and Sterman (see Refs. mm and references therein). Recently the advent of the soft- 
collinear effective theory (SCET) [3H5] facilitates the proof in various scattering processes 
with the systematic treatment of power corrections. For example, the factorized dijet cross 
section in annihilation can be schematically given by the convolution of the hard, 

collinear and soft parts as 


a = ® Jnifi) ^ ( 1 ) 

where H is the hard coefficient, and Jn are the collinear jet functions which describe the 
final-state energetic collinear particles in the lightlike directions n and n. The soft function 
S describes the soft interaction. In pp scattering, the scattering cross section is extended to 
be given as 

cr = H{Q^,p)Jn{p) (8) Jn(/i) 0 fp{xi,p) (g) fp{x 2 ,p) 0 5'(/i), (2) 

where the additional fp comes from the parton distribution functions for initial hadrons. 

Once the factorization theorem is established, each part can be computed perturbatively 
and can be resummed by solving the renormalization group equation. However, there are 
many constraints for the theoretical computation to be consistent in order to satisfy the 
factorization theorem. First, each part should be infrared safe. Otherwise, the factorization 
itself does not have any physical meaning. Secondly, since the scattering cross section is 
independent of an arbitrary renormalization scale, the sum of the anomalous dimensions of 
all the factorized parts should cancel at a given order in perturbation theory. 

In high-energy processes, final-state hadrons usually form a collimated beam which is 
called a jet. Jets provide important information on the physics within and beyond Standard 
Model. In describing high-energy collisions in terms of jets, the factorization property is also 
essential to theoretical predictions. There are diverse ways to define jets depending on the 
kinematics of the scattering, the detector design, etc.. The jet definitions are represented 
by jet algorithms such as the Sterman-Weinberg (SW) algorithm |6], the cone algorithm [7], 
the JADE algorithm [S], the algorithm [H], the anti-ZcT’ algorithm[Tn], to name a few. 

The first theoretical attempt to implement the jet algorithm at the parton level is the 
Sterman-Weinberg jet. The issue is how to combine adjacent particles into a jet without 
ambiguity. The important ingredients for a good jet algorithm are that it can be employed 
conveniently both in experiment and in theory. And it should be infrared safe. However, 
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implementing a convenient jet algorithm in experiment is sometimes difficult to realize in 
theory, and vice versa. 

In this paper, we hrst establish the factorization theorem in the dijet production from 
electron-positron annihilation in the framework of SCET. The cone-type algorithm suggested 
by Ellis et ah m and the Sterman-Weinberg algorithm are employed as specihc examples. 
We explicitly show that each factorized part is IR safe at order ag. This is verihed in two 
ways. First we compute each factorized part, using the dimensional regularization both for 
the ultraviolet (UV) and the infrared (IR) divergences. And we can put a rapidity regulator 
in the collinear Wilson line to regulate the IR divergence, while the UV divergence is still 
handled with the dimensional regularization. 

We consider two kinds of jet functions. One is the unintegrated jet function, and the 
other is the integrated jet function which is obtained by integrating the unintegrated jet 
function with respect to the invariant mass squared of the collinear part. Intuitively, the 
integrated jet function should be obtained by integrating the unintegrated jet function, 
which is true in our calculation. Related to this topic, different kinds of jet functions have 
been discussed in Ref. m- The measured jet function describes the jet with thrust, while 
the unmeasured jet function describes the jet without measuring the thrust. However, the 
unmeasured jet function is not obtained by simply integrating the measured jet function. 
This may look confusing at hrst sight, but this results from different ways of power counting 
in the measured and unmeasured jet functions. And in deriving physical results such as the 
scattering cross section, using the jet functions in both methods are consistent. This will be 
explained in detail. 

The idea of applying SCET to the factorization of the dijet cross section has been tested 
in various directions. An interesting attempt was initiated by Cheung et ah [12] where the 
authors carefully analyzed the phase space according to jet algorithms. The SW, JADE, and 
kr jet algorithms were analyzed and the dijet cross sections were computed. Though they 
did not explicitly state the factorization theorem, they implicitly worked in the framework 
of SCET by computing the collinear and the soft parts. 

Jouttenus PI focused on the jet functions in e’''e annihilation with the factorization 
scheme in SCET with the SW jet algorithm. In Ref. PI, the jet cross section was computed 
by deriving jet operators and computing the matrix elements of those operators. In Ref. HH, 
jet shapes were comprehensively analyzed employing SCET. They used a cone-type jet 
algorithm and a kr algorithm as examples to explore the properties of the “unmeasured” 
(corresponding to “integrated” in our terminology) and the measured (unintegrated) jet 
functions, as well as the soft functions. They pointed out that the unmeasured jet function 
is not obtained by integrating the measured jet function because a different power counting 
was involved. This will be discussed in detail compared to the fact that our integrated jet 
function is obtained by integrating the unintegrated jet functions. 

The structure of the paper is the following: In Section 2, a factorization theorem for 
the dijet cross section in annihilation is presented. In Section 3, the characteristics 

of the generic cone-type and the SW jet algorithms are described, and the available phase 
space for the collinear and soft parts is clarihed. We also delineate the power counting of 
the various quantities appearing in the factorization theorem. In Section 4, the integrated 
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jet functions with both jet algorithms are computed using the dimensional regularization 
both for the UV and IR divergences, and the unintegrated jet functions are presented in 
Section 5. In Section 6, the soft function is computed using the dimensional regularization 
to order in both jet algorithms. In Section 7, we verify the SW dijet cross section in 
SCET is the same as that in full QCD at order and show that the anomalous dimensions 
of the factorized parts are summed to cancel. In Section 8, we discuss the relation between 
our integrated, unintegrated jet functions and the unmeasured, measured jet functions and 
explain why both approaches are consistent. Finally in Section 9, we give a conclusions and 
perspectives. In Appendix, the jet function in the SW jet algorithm is computed using the 
rapidity regulator, in which the UV divergence is handled by the dimensional regularization, 
while the rapidity divergence is handled by the rapidity regulator. 


II. FACTORIZATION OF THE DUET CROSS SECTION IN SCET 


In annihilation, as well as in hadron-hadron scattering, a great number of jets may 
be produced at high energy. A jet algorithm should be implemented to identify a scattering 
event with N jets. A proper jet algorithm should separate each collimated beam of particles, 
and merge nearby jets without ambiguity with infrared safety. However, in this paper, we are 
interested in the conhguration in which there are two back-to-back jets with soft particles, 
which is called a dijet event. The property of V-jet events, for example, N jettiness on event 
shape, can be found elsewhere 

We consider the dijet cross section in —)■ jij 2 X, where ji and j 2 describe collimated 

beams of particles. Because ji and fy are almost back-to-back in the center-of-mass (CM) 
frame, we describe collinear particles inside ji and fy in terms of the lightcone vectors n and 
n respectively with = 0, n -n = 2. The dijet scattering cross section is given by 

cT(e+e“ -^jij 2 X) = - PxJ(0| |X) (X| J^|0), (3) 


where is the center-of-mass squared of the e’''e pair. Here ao is the Born cross 

section for a given flavor / of the quark-antiquark pair, given by 


ao = 


At^o^QINc 

3Q^ 


( 4 ) 


The hnal states denoted by X consist of two groups of collinear particles in the n and n 
directions respectively, and the remaining soft particles. 


|A) = |X„)^|X^)^|A,). (5) 

The factorization theorem can be proved in terms of the hadronic state |X) without this 
decomposition [16] , but we use Eq. (|^ for simplicity. The collinear and soft particles will be 
combined into jets according to the jet algorithms, which are described in the next section. 
In SCET, the collinear particles in the n and n directions have the momentum scaling as 

Pn = in-pn,Pn,n-pn) = {p~,pi,Pn) ~ <5(1, A, A^), 

~<5(A^A,1), (6) 
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where A is a small parameter in SCET, while the soft momentum scales as 


= {Ps .p^.pt) ~ A^ A^). 


( 7 ) 


For the photon exchange, the electromagnetic current is written at leading order in SCET 
as 

= ( 8 ) 

where Xn is a gauge-invariant collinear quark held combined with a collinear Wilson line, 
(x) is the soft Wilson line dehned as HU 



P exp 


W 



dsn ■ As{sn) 


(9) 


where ‘P’ denotes the path ordering. Note that the specihed path in W is different from the 
standard case of W(x) |5], where the path is given by [—oo, x] along the lightcone vector n. In 
Eq. (|^, C{Q‘^, fi) is the Wilson coefficient obtained from the matching between full QCD and 
SCET. The hard coefficient in the scattering cross section is given by /i) = |C(Q^,/i)p, 

which is given to one loop as [T8] 


H{Q\p) 


1 + 


(^sCp ( . 2 




( 10 ) 


We dehne the (unintegrated) jet function in the n direction as 


|:(0|x:|V„>{Wlx^|0) = / ^n-Px,.|x(pL.A)V'’, 


( 11 ) 


while the jet function in the n direction is obtained by switching n and n. The jet function 
at lowest order is normalized as Jn\p^) = ^{p^)- Because the collinear fields do not interact 
with the soft helds any more after the held redehnition |5], the scattering cross section can 
be written as 

a = aoH J d'^pid'^p2^^Jn{p\,p)Jn{pl,p) 

X 1 Tr(0|f4y;|W)(x,|yiw|0), ( 12 ) 

where we put pi = px„, P 2 = and ps = pxs- AVe choose the coordinate system such 
that Pi = 0, then the phase space integral for pi can be written as 


d'^pi = 2 tt J dpi dpi 


-^^+(Pi -pT 
2 


TT 


dPidpiipi 


-^2 


(13) 


In the CM frame, the delta function in Eq. (|12|) is rewritten as 
S{q-Pi-P 2 - Ps) = 2S{Q-pi - P 2 


Ps )d{Q - pt - pt - pt)d^^\pi: + P 2 + pi)- (14) 


Applying the power counting in Eq. (|^ we can simplify further the delta function in Eq. (14). 
If we consider a differential scattering cross section for observables of order A^, the soft part is 
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correlated with the collinear parts. For example, if we look at some distribution concerning 
the range of 0{\^) in the p~ component, we have to keep P 2 rsj 0{X^) in 5{Q-p^ -P 2 -p ,) 
when Q — Pi ^ 

Unless we consider the pt distributions of O(A^), we can suppress in + + 

since p ^2 ~ ^(-^)- Then the scattering cross section becomes 

a = aoH{Q^) J dp'tdpidptdp2 ^^^^^‘^ Jn{pl)Jn{pl) (15) 

X ^5(g - Ph - P2 - p:)S{Q -Pt-Pt- P+)^Tr(o|y^y;|x,)(x,|yTyk|o), 

= aoH{Q^) j dpfdpidp^dp^dl'^dr^^^^^^^diQ - pf - P 2 - 1~^)6{Q - P 2 - Pi - l~) 

X J„(p?) J^(p^) E ^Tr(0|y^y;|x,)(x#(/++ iP+)5(r + iP7)y;tyk|0), 

= aoH{Q^,p) J dp2dp'lp];ptJn{pl,p)Jn{pl,p) J dl^drs{l^,r, p), 

where pf = Q — pt — P 2 = Q ~ Pi ~ ^ Pi = pt — function S'(/+, l~) is 

dehned as 

S(i\ n = iTi'{0|y^pi(;+ + VPHI- + P;)h'n|o>. (le) 

If a jet algorithm is applied, it constrains both the collinear and soft parts. The total 
dijet scattering cross section can be written as 

o-j = c^oi^(g^/i) J dptdp2Q‘^Jn,e{pl)Jn,0{pl) j dtdrSe{l^,l~) 

= croi/(g^/i) J dplJn,@{pl,p) J dplJn,e{pt p)<Seip) 

= o'oH{Q‘^, p)Jn,e{p)>Jn,e{p)<Se{p), (17) 


where the subscript 0 means the jet algorithm is implemented. And Sq is the integrated 
soft function given by 

Se{p) = J dl+dl-Se{l^,l~,p). (18) 

If we do not apply jet algorithms, the integrated soft function becomes 1 to all orders in 
as because of the unitarity. For later use, we call Jn,e{p‘^,p) the unintegrated jet function, 
and Jn,@{.p) = / dp^ Jn,@{.P^^ p) the integrated jet function. When we focus on the dijet total 


scattering cross sections without specihc observables, the factorization theorem in Eq. (17) 
is enough. 


As clearly seen in Eq. (17), the scattering cross section is factorized into the hard, collinear 
and soft parts. And since the jet cross section is a physical observable, it should be inde¬ 
pendent of the renormalization scale p. Also the double differential scattering cross section 

daj 


2j 2 = cro77(g , p)Jn,e{Pa P)Jn,eiP 2 ^ P)'Seip) 


dpidp\ 


(19) 


should be also independent of p. That is, the integrated jet function and the unintegrated 
jet function have the same anomalous dimensions. 
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Note that the double differential scattering cross section is not physical since pf and pi 
are not physical observables, but it can be used as an intermediate step to obtain the total 


jet cross section. However, pf, pi in Eq. (19) are not the measured invariant jet mass squared 


with the corresponding jet algorithm. The physical invariant jet mass comes from the sum 
of the momenta of the collinear and soft particles inside a jet. The invariant masses squared 
for the n and n jets when they contain a soft gluon are given as 


— [pi + lY — p\ + QI+ + Cl(A^), 

= (p2 + 0^ = P 2 + Q^- + 


( 20 ) 


And the double differential scattering cross section with respect to these physical invariant 


jet masses is obtained from Eq. (15) as 
daj 




= croi^(Q^/i) dpldpldl+dl Jnpipl, p)Jn,e{pl, p)S0{l^,l ) 


( 21 ) 


X 6 -pI- - pI - Ql~ 

= croi/(Q^/i) J - (5/+, p)7^,0(771^2 - Q/“,p)S'e(/+,/“,/!). 


This factorization theorem is similar to the result in Ref. [19], where the double differential 
scattering cross section for top jet mass has been studied. 

And the single differential scattering cross section is given as 

da 


drnj^ 


Y = (ToH{Q‘^,pt) j dpjdpldl^dl Jn,e{pl, p)Jn,eipl, J - pf - Ql^'^ 

= o■oi^(Q^/u) J - Q/+,/i)S'e(/+,p) 54 ,e(p), ( 22 ) 


where S'©(/■*■,/r) is given by 




(23) 


If we are interested in the dependence of the invariant jet masses, Eq. (21) or (22) should 


be used. They are given by the convolution of the jet and the soft functions. But if we are 
interested in the total jet cross section, Eqs. 0 and ( p)^ can also be used. In order to 
obtain the integrated and unintegrated jet functions, Eqs. (17) and (19) are employed from 


now on. Note that the total dijet cross section in Eq. ( [ItI ) is a product of of the integrated 
jet functions and the soft function. 


III. JET ALGORITHMS 

Among many jet algorithms, here we describe a generic cone-type and the SW jet algo¬ 
rithms. The characteristics of other jet algorithms will be treated elsewhere. The basic idea 
of a jet algorithm is to combine final-state particles into a jet if those particles are within the 
radius of the jet cone R (cone-type algorithm) or the cone of half-angle 5 (SW algorithm). 
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FIG. 1. A schematic diagram for the jet function at next-to-leading order. The jet momentum is 
p, and particles in the loop carry momentum I (gluon) and p — I (quark). 

The main difference lies in the fact that the cone-type algorithm is concerned about the 
angle of a particle with respect to a specihc axis such as the thrust axis of the jet, while the 
SW algorithm is concerned about the relative angle between the particles to combine them 
into a jet. Since the algorithms are similar conceptually in the sense that particles form a 
jet inside a jet cone, they show similar behavior in the jet and the soft functions. 

Let us be more quantitative on the jet algorithm. We consider at the next-to-leading 
order, and there are two particles at most in a jet. In e’''e“ annihilation, a qq pair will be 
produced at leading order forming a back-to-back jet. A gluon can be emitted at order a<j, 
and it can be combined with a quark or an antiquark to form a jet. A qq pair can also form a 
jet, but its contribution to the jet cross section is suppressed compared to the former cases. 

We show schematically the assignment of momenta in Fig. [Tj The total momentum of the 
hnal state is p with p_ = Q for the jet in the n direction. And p^ is the invariant-mass squared 
of the collinear jet. In order to obtain the scattering cross section, we cut the diagram in all 
possible ways and add their contributions. When we cut the loop, there are two hnal-state 
particles with the momenta I (for a gluon) and p — I (for a quark or an antiquark). If we 
cut a single line, it gives virtual corrections. The collinear gluon momentum in the n 
direction scales as 

P = (/_,/x,4)~Q(l,A,A2), (24) 

where Q = P- is the largest momentum component of the collinear jet. 

A cone-type algorithm collects all the particles within a fixed cone of radius i?, which 
determines the jet size. The jet axis is determined, for example, by a thrust axis. There 
are many issues on how to construct cone jets including infrared safety. But if we conhne 
ourselves to the jets at the next-to-leading order, there are at most two particles in a jet 
and all the cone-type algorithms are reduced to a generic cone-type algorithm. For collinear 
particles, if two particles are in the cone of radius R, that is, if Oq < R, and 9g < R, where 
6q and 9g are the angles of the quark and the gluon with respect to the jet axis, they form 
a jet. In terms of the momenta shown in Fig. this corresponds to the cone algorithm for 
the collinear particles: 

©cone = 0 0 , (25) 

8 







where t = tani?/2 with the jet cone radius R. Each term in ©cone represents 9g < R and 
9g < R respectively. 

In the SW jet algorithm, a two-jet event is dehned as the one in which the total energy Q 
except the fraction (3 is deposited in two jet cones with half angle 5. This algorithm can be 
split into a collinear part and a soft part, according to which the factorized result is written 
in Eq. For the collinear part, the two particles are energetic and, if we let 9 be the 

angle between the two particles, the jet algorithm states that 9 < 26, which is equivalent 
to 1 — cos 6* < 1 — cos 25 for 0 < 9 < tt. For the jet in the n-direction, combining the 
on-shll conditions R = 0, {p — = 0, and using the hierarchy of lightcone momenta, this 
corresponds to 

Q% 


©SW — > 




(26) 


L(Q-L)2. 

with t = 6. 

We will use t both for the cone radius {t = tani?/2), and the angle {t = S). So far, the 
cone size is independent of the power counting of momenta in SOFT. For example, if R = tt 
for the entire hemisphere, t = tani?/2 ~ 0{1). However, if t is too small, the constraint 
of the jet algorithm does not work. The size of the parameters t and f3, related to the jet 


algorithm, should be determined by the physics we consider. As can be seen in Eq. (26), 
with the fact that the collinear momentum scales as (Z_,Z_l,Z+) ~ Q(1,A, A^), the power 
counting R ~ A^ produces a nontrivial jet algorithm. Therefore here we consider the case 
with f ~ A, and t can be used interchangeably with tani?/2 or 5. For the jet in the n 
direction, the corresponding condition is given by switching Z_|_ and Z_. 

To avoid double counting, the zero-bin contribution should be considered m, where the 
collinear momentum R approaches the soft-momentum limit 


/^ = (Z+,Zx,Z_)~g(A^A^A2). 
The corresponding zero-bin phase space is expressed as 

L 


e(“) = eif > 


L 


(27) 


(28) 




FIG. 2. Constraints on the phase space of the collinear jet function. The shaded regions are allowed 
by the jet algorithms, (a) Phase space in the cone-type jet algorithm, (b) Phase space in the SW 
jet algorithm, (c) Zero-bin region for both jet algorithms. 
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FIG. 3. Phase space available for the soft function. 


This is common to the cone-type and the SW jet algorithms. The phase spaces given by 
the jet algorithms are shown in Fig. Fig. I (a) is the phase space available for the naive 
collinear contribution in the cone-type jet algorithm, Fig. (b) is the phase space available 
for the SW jet algorithm, and Fig. (c) is the phase space for the zero-bin contribution for 
both jet algorithms. 

We use the same constraints for soft particles in both algorithms. That is, when a soft 
gluon is inside the cone either in the n or n directions, it is included in the jet. An additional 
constraint on the soft function is about a jet veto. If a soft particle is outside the jet, the 
energy of the soft gluon is less than (3. The soft momentum scales as QX^, and the nontrivial 
constraint from the jet veto occurs when (3 ~ A^. From now on, we £x the power counting 
on the parameters in the jet algorithm as t ~ A and ~ A^. 

Then the constraint on the soft function from the jet algorithm can be expressed as 


©soft - 


©(/+ < 2(3Q), (outside the jet). 


(29) 


Fig. shows the phase space for the soft function. The shaded area is allowed by the jet 
algorithm in Eq. (29) 

Note that we include the constraint from the jet veto in the jet algorithm of the soft 
function, not in the collinear function. This is obvious from the power counting since /? ~ A^ 
and 1+ + I- ~ Q for the collinear gluon, hence there is no veto needed for the collinear 
particles. Instead the zero-bin subtraction eliminates the soft limit in the collinear jet 
function. However, in Ref. 1131 . the author considered only the jet function with the same 
power counting on t and /?, but this jet veto was also included in the jet function. This 
seems to be a different dehnition of the jet function, compared to ours. But the jet function 
in Ref. [13] does not depend on (3 after the zero-bin subtraction. That is, the inclusion of a 
jet veto in the jet function does not alter the result. However, we stress that our definition 
of the jet algorithm strictly respects the specified power counting of the parameters. 

A question arises if the cone-type and the SW jet algorithms can also be employed at 
hadron colliders. The main difference between e’''e“ annihilation and pp scattering is that 
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the center of energy is fixed in e'''e“ annihilation, while it is not hxed in pp scattering. And 
the center-of-energy frames at the partonic level in pp scattering have a boost invariance 
along the beam direction. Therefore the jet algorithm in pp scattering should be dehned in 
a boost-invariant way, for example, in terms of the rapidity and the azimuthal angle. The 
cone-type algorithm can be employed in pp scattering, and the method of dehning the jet 
axis should be employed. The SW jet algorithm seems to be specihc to e’''e“ annihilation. 
However, at the next-to-leading order in which there are three hnal-state particles, the 
inclusive k^p algorithms {kp, Cambridge/Aachen, and anti-fc^’) reduce to the same constraint 
as the SW jet algorithm, and can be used in hadron collisions. 


IV. INTEGRATED JET FUNCTION 


The unintegrated jet function in the n direction is dehned as 

II(O|Xnl^n)0j(^n|An|O) = J 

where 0j is the jet algorithm. The integrated jet function Jn,e{.k) is dehned as 


(30) 


Jn.a(li) = J dp^Jn,e(P^< f*). 


(31) 


where p^ is the invariant mass squared of the collinear jet. The unintegrated jet function is 
the collinear jet function with p^ hxed, and it is expressed in terms of distribution functions. 


When integrated over p^, the integrated jet function is obtained. Eq. (31) should hold to all 


orders in as- But in current literature, there is some confusion between the integrated jet 
function, which is obtained by integrating the unintegrated jet function, and the inclusive 
jet function in which there is no jet algorithm involved. 

In a diherent context in Ref. HU, the unmeasured jet function (which is the same as the 
integrated jet function in our case) is not obtained from the measured jet function (which 
corresponds to the unintegrated jet function, but diherent) by integrating with respect to 
the thrust. The authors claim that the small parameter in power counting is diherent. That 
is, in a measured jet sector the small parameter is ~ A, while R ~ A°. Indeed diherent 
power counting yields diherent answers, but as far as the total jet cross section is concerned, 
the result through the paths from the integrated jet functions, and from the measured jet 
functions is the same. The relation between these two approaches and the physics behind 
them will be clarihed in Section VIII. 

In Sections IV and V, we explicitly compute the integrated and the unintegrated jet 
functions separately. The diherence lies in the diherent order of integration of the delta 
functions, but there are subtleties which arise in the detailed calculations. What we also 
focus on is whether the factorized parts, the jet function and the soft function, are infrared 
hnite themselves. If they contain infrared divergences though the sum of the jet and the 
soft functions do not, the factorized result is not physically meaningful. In order to show 
the IR safety explicitly, we employ the dimensional regularization with D = 4 — 2e and the 
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FIG. 4. Feynman diagrams for the jet functions at one loop, (a) virtual correction, (b-c) real 
gluon emission. The mirror images of (a) and (b) are omitted. 


MS scheme, and regulate both the UV and the IR divergences. In the MS scheme, we put 
4vr/i|g =/iV^. 

We can also employ the dimensional regularization for the UV divergence and the rapidity 
regulator only in the collinear Wilson line [211 122 ] to regulate the rapidity divergence as 


1U„ = ^ exp 


9 


perm 


n ■ P + A + fO 


n ■ Ar^ 


(32) 


The hnal collinear results with the zero-bin subtraction are independent of the regulator 
A though it appears in individual diagrams, and the result is the same as that with pure 
dimensional regularization. We present the result for the integrated jet function in the SW 
jet algorithm in Appendix for completeness. 

The Feynman diagrams for the jet function, say, in the n direction, are shown in Fig. 
The dashed line represents the cut, therefore Fig. |^(a) is the virtual correction, while Fig. 
(b) and (c) are the contributions from the real gluon emission. The mirror images of (a) and 
(b) are omitted in Fig. but they are included in the computation. And we also include the 
wave function renormalization and the residue from the self-energy diagram for the quark. 
Each diagram is accompanied by the corresponding zero-bin contribution, which should be 
subtracted to obtain the collinear jet function. 


A. Cone-type jet algorithm 


The naive collinear virtual correction in Fig. (a) is given by 




d^l Q-l- 

(27r)^ /2(p -/)2/_ 


a.C 


k fiTi ' ^ t 


2n ^euv eiR^^eiR 


-M -h In 


/i 




(33) 


where Q = p_ is the largest component of the lightcone momentum in the n direction. In 
obtaining the hnal result, we use the fact that 




duu 


-l-e 


^UV 


1 

Gr’ 


(34) 


where u is the integration variable with the dimension of mass. In the dimensional regular¬ 
ization in which the UV and the IR poles are not distinguished, the above integral is simply 
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zero because it is a scaleless integral. However, if we regulate the IR divergence by inserting 
a regulator as 


j/ / du{u + 5 ) 


-l-e 


<;uv 


( 9 ' 


<;uv 


1 ^ 


(36) 


the UV pole can be explicitly extracted. Therefore in order for the integral in Eq. (34) 
without the regulator to satisfy the requirement that it be zero, the pole structure should 
be of the form given in Eq. (34). Note that we presume euv > 0 to assert that the integrand 
at inhnity vanishes in Eq. (35). For the same reason, we presume that cir < 0. 

The corresponding zero-bin contribution is given as 


= 


a.^C ^ ^ ^ 


, -)t 

2% ^euv 


The collinear contribution from Fig. (a) is given by 
Ma = Ma 


, QaC F / 1 

M ° ( 


27r ^euv ciR^^euv 





(36) 


(37) 


Since the virtual corrections are independent of the jet algorithm at order Og, this result is 
also true in the SW algorithm to be described below. 

The naive collinear real gluon emission in Fig. (b) is affected by the jet algorithm and 
is given as 


Mh = 




27r r(l — e) 


Q 


(38) 




27r r(l — e) 




fQ Q-l 1 r{Q-l-)hyi. 

Iq/2 Q Jo 


agCp r 1 , 1 , 1 1 




27r L2e2j^ ejR V '2 QH^ 4 QH^ 


+ +bn=4hT + 2 + 21n2-|-7r'^ 
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As can be seen clearly in the integral, there are only IR poles in this naive collinear contri¬ 
bution. 

And the corresponding zero-bin contribution is given as 


0 _ oisCp 




27r r(l — e) 

271 r(l -e)Jo - Jo + + 


rt^L 


(39) 


Care must be taken in computing this integral in dimensional regularization. If we naively 
perform the /+ integral, it yields —with the presumption that e < 0 to guarantee 
that the integral vanishes at /+ = 0. That is, the divergence is of the IR origin. This 
prescription is hne for the /_ integral near /_ ~ 0. However, when we perform the /_ 
integral, the integral should also be regulated as /_ approaches inhnity. It means e > 0 in 
this integral and the presumption that e < 0 in the integral is violated as L_ approaches 
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infinity. Therefore we cannot obtain a meaningful answer if we naively compute the integral 
as it is. 

The trick to avoid this problem is to write the integral as 


dl.l 


-l-e 






di.r^-^+ I duz 


-l-e 


= / dlJ 


-l-e 




di+iz^-^+ I dUZ 


-l-e 


dlA 


-l-e 






-l-e 


dlA 




(40) 


where r] is some positive energy. In the hrst line, the hrst integral on the right-hand side 
yields the IR poles (e < 0). The second integral, as it is, is inconsistent because the /+ 
integral forces e < 0, while the integral on l_ requires e > 0 as /_ —>■ cx). Therefore the 
second /+ integral is decomposed to yield the final result. The last integral is purely of the 
UV origin, while the second integral has the UV origin from the /_ integral, and the /+ 
integral can be treated by Eq. (34). With this careful manipulation of separating the UV 
and IR regions, the hnal result is independent of the arbitrary scale 17 , and the zero-bin 
contribution is given by 




agCp 

2ti 


[9 

^ 1 

1 x2 1 

-+ 0 

^ 1 

)lnf2l 

L2 

Vuv 

^IR^ 2 

Vuv 



(41) 


If this trick were not used and if the integral were computed carelessly, the result would 
be totally different and we get a meaningless jet function. The collinear contribution from 
Fig. (b) is finally given by 


M, = M,- = 


0 _ C^sCf 


271 


1 1 1 . /i\ 1 


+ 




2euv cuvUr ^ir ^ cuv 


Inf 


In ^ + 1 ln 2 ^ + 2 + 2 In 2 - -tt^ 

24 


Q2t2 4 Q2t2 

The naive collinear contribution shown in Fig. (c) is given 


(42) 


as 




27r 

agCp 

271 


' dl+dlJ^l-^-^Q, 


r(l-e) 

r /■Q /2 


(43) 


^^T{l-e)Uo 

( _ ^ _ 

47r V ejj, ^Q2t2 


n2 I 
Jo 

1 - 21n2V 


1 /-Q I 

dl+lZ^-^+ / dl- 


l-e 




Q/2 Q2 Jq 


dlA 


-l-e 


Note also that the pole is of the IR origin. The zero-bin contribution is suppressed at this 
order. 

Finally, the wave function renormalization and the residue from the self-energy correction 
of the fermion at order ag are given as 




C^sCp / 1 
dvr Vuv 



(44) 
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Combining all the contributions by including proper mirror images, the total contribution 
to the integrated jet function is given by 


d^coll — 2(iffa + Mb) + Me -\- 


hi) 


(45) 


r 1 


271 


+ 


1 /3 


^e 


uv 


euv 




2 ^ QH^ 




^ 2 QH^ 


. 1,2 

^2 QH^ 


7 Stt^ 

H-h 3 In 2- 

2 12 


Note that the integrated jet function in the cone-type algorithm is indeed IR hnite. By 
removing the UV poles, the integrated jet function in the cone-type jet algorithm at order 
as is given by 




I 

QH^ 2 


^ 7 

3 In 2 


12 


(46) 


This coincides with the unmeasured quark jet function in Ref. m- 


B. SW algorithm 


Since the SW algorithm and the cone-type algorithm are essentially the same except the 
fact that the cone is determined either by the relative angle of the two partons, or by the 
angle with respect to a hxed jet axis such as the thrust axis, the calculation for the jet 
function in the SW jet algorithm is very similar to the jet function in the cone-type jet 
algorithm except some constants. 

Since the virtual corrections in Figs. (a) and the wave function renormalization are 
independent of the jet algorithms, they are the same as in the case of the cone-type jet 
algorithm. The naive collinear real gluon emission in Fig. (b) is affected by the jet algorithm 
and is given as 


Mb = 


asCp (ei'"/r^)^ 
271 F(1 — e) 


asCp (eT'E/r^)^ rQ Q- l_ ri-t\Q-i.)VQ' 


277 F (1 — e) Jo 


dl- 


Q 


dij 


-l-e 


agCp 




H- (1 + X In 




+ In 


JJ— + 1 In2 JJ— 

QH^ 4 QH^ 


2 ,32 

-t- 4 — -TT^ 


(47) 


As can be seen clearly in the integral, there are only IR poles. Since the zero-bin constraint 
is the same as that of the cone-type jet algorithm, the zero-bin contribution is also given by 
Eq. (41). Finally, the collinear contribution from the real gluon emission is given by 

OtsCp 


Mb = Mb- MJ = 


277 


1 1 
+ 


9 2 ' + -+ 

2euv euv^iR eiR ^ Q 


Int 


euv 


In 


+ i ln2 

QH^ 4 QH^ 


-|- 4- 77 
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(48) 
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The naive collinear contribution shown in Fig. (c) is given as 


;l-e 


Mr = ^^(1 - [ dl+dlJ-^l-^-^Qsw 


2n 

agCp 


r(i-6) 




(49) 


( 1 -^) 


(e^'E/i^)" fQ Iz^-^ A-t^Q-i.)yQ^ 


2tt '^r(l-e)7o Q2 

sCf ( _ 1 _, _ Q 

47r V eiR QH^ 


dlS^^ I dlJZ^-^ 


Note also that the pole is of the IR origin. The zero-bin contribution is also suppressed at 
this order. 

Adding all the contributions including mirror images of (a) and (b), the total contribution 
to the integrated jet function for the SW jet algorithm is given as 


Afcoii — 2,[Ma Mb) Me -\- -\- R'^ 


?(i) 


(50) 


exqC F r 1 
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1 /3 


(; 


-|- In 




3 /i^ 

^ 2 ^ 


1 1 2 

^2 ^ 


13 3 . 

+ w - tTT^ 


27r Le^v ' <^uv'"2 ' 

Therefore the integrated jet function in the SW jet algorithm at order is given by 


,/^)- ^ (9 ’^02+2 + 9 02+2+ 9 4^)' 


27r 


QH‘^ 




(51) 


This is also consistent with the /cr-type result in Ref. m, since the SW algorithm and the 
fc^-type algorithm are the same at this order. Note that, compared to the contribution to the 
integrated jet function in the cone-type jet algorithm, all the divergence structure and the 
logarithmic terms are the same. The only difference comes from the constant terms. This 
is obvious when we look at the phase space for both jet algorithms Fig. As /_ goes either 
to zero or Q, the phase space coincides, which yields the same divergences and logarithms. 
The difference of the remaining phase space yields different constant terms. That is why we 
have stated that the cone-type and the SW jet algorithms are similar. 


V. UNINTEGRATED JET FUNCTION 

Here we derive the unintegrated jet function for the cone-type and the SW jet algorithms. 
It should yield the integrated jet function when integrated over In computing the in¬ 
tegrated jet function, we integrated the delta function responsible for the on-shellness of 
the final-state particles with respect to hrst, and the remaining and /_ integrals were 
performed. Here we do not integrate over Instead, we perform the /+ and /_ integrals 
hrst and express the result in terms of the distribution functions with respect to p^. The 
integrated jet function should be independent of the order of integration, and our calcula¬ 
tion explicitly shows that fact. However, there are some caveats on how to deal with the 
distribution function, which will be explained here. 
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A. Unintegrated jet fnnction in the cone-type algorithm 


The Feynman diagrams for the nnintegrated jet fnnction are also given by Fig. The 
virtnal correction in Fig. (a) and its zero-bin contribution are the same as the result in 
the integrated jet function except the delta function and the collinear result is given as 


Ma = Ma 




( ^ 


( ^ 

+ In ^ + l] 


Ak 


Q ^ 


(52) 


The naive collinear contribution in Fig. (b) is given 

ttsCp 1 

271 F(1 — e) 


as 


Mb = 


dl+dl-6{{l - ^ ^~0 c 




277 F(1 — e) (p^ji+c 


dyy '(1 - y) 


l-e 


(53) 


where y = l-jQ and + p^. Note that p^ < Q‘^t^ in order for this integral to 

exist. Since p^ also appears in the integration limits, the dependence on p^ does not come 
only from in front of the integral. However we observe that Eq. (53) is singular at 

p^ = 0, and can be regarded as a distribution function. It can be written as 


Mb = A6{p‘^) + [H(p2)]M2, 


(54) 


where all the singularities are concentrated at p^ = 0 and the remaining part gives a distri¬ 
bution function away from p^ = 0, which we call the M^-distribution. It is defined as 




dpW)Wfip")= / dp^g{p^)f{p^)- / dp^p(p^)/(0), (55) 

Jo Jo 


for a well-behaved test function /(p^). Here p(p^) is a function which diverges at p^ = 0 and 
A is an appropriate finite upper limit for the relevant physical processes in consideration. It 
has the following properties: 


rM^ 


dp^[g{p^)\M^f{p^) = / dp‘^g{p'^)(^f{p^) - /(O)), 


rM^ 


dp^[g{p‘^)]M2 = 0 . 

For the jet algorithms here, we can put M = Qt in the final result since p^ < Q^t^. 
The M^-distribution part can be computed with e = 0, and is given as 


(56) 


rQ'^PluP- X 

dy- 


y 


— p^ 

+ p 2 


In 




pZ 


p^ Jp'^/w^ " y p^ 

The part proportional to 5(p^) is obtained by integrating Mb with respect to p^ as 

rM^ dp2 fQ^P/w^ 


(57) 


/o {p‘^y+^ 


dyy "{i-yY 


(58) 
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Using integration by parts, the calculation is straightforward, and the hnal result is given as 


A/,= ^P(P) 


271 


^ 5 —I-(1 H— In — 

6irV 2 


, + -ln^^ 






In —--In 

M 2 2 


I r) 2+2 




2+2 


- 21n 




.U(- 


M2 QH'^ + M^i Lp2V QH^+p'i 

The zero-bin contribution is given by 


Q2t2_p2 q2^2 

+ In —^ 


M 2 


}.(59) 


MS = 


agCp 


277 r(l - e) (p2)i+e yp 2 /Qt 2 Q \QJ 


v/p 2 / Qfi Q ^ 


dl- fl- \ - 1 -e 


— 


(60) 


We have to be careful in expressing this as a distribution function. Note hrst that the 
upper limit of the integration is set to inhnity. In the naive collinear contribution, the 
collinear momentum scales as (/_, U, /+) ~ Q(l, A, A^). In the zero-bin limit, scales as 
(/_,/_L,/+) ~ Q(A^,A^,A^). The upper limit is of order 1 since it lies in the range 1/2 < 
Q‘^t^/+ P^) < 1 for 0 < while 1-/Q ~ 0{X^). Therefore in the effective 

theory, the upper limit is set to inhnity. It also means that p^jQ’^ is at most 0{\^) or 
smaller. If we assume that p^ is of order as in the naive collinear contribution, the 

lower limit of the integral in Eq. (53), p‘^/{Q‘^t^ X-p"^), becomes of order 1 since ~ (!1(A^), 


while 1-/Q ~ 0{X^). In this case, the lower limit should be set to inhnity also, which is 
not sensible. Therefore, in computing the zero-bin contribution, we should consider p^ much 
smaller than Q^X^, say, Q^A^. It means that the zero-bin contribution is concentrated near 
the origin, and there is no remaining distribution function away from = 0 which describes 
the region p^ ~ Q^X^. In summary, the zero-bin contribution is proportional only to 5(p^). 

With the above argument, the coefficient of S{p‘^) in the zero-bin contribution can be 


obtained by integrating Eq. (60) with respect to p^ from 0 to inhnity. Here also, the trick 


similar to Eq. (40) should be employed to avoid using the wrong sign of e. By shuffiing the 


integration region appropriately and introducing an intermediate scale rj, we can obtain the 
unambiguous UV and IR divergent terms as 


M° = 


agCp 

~2YY 
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2 1 
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1 
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+ 

>1 

2 1 

1 

1 

c-f* 

to^ 
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'euv 

Ur^ 2 

'euv 
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(61) 


Note that the coefficient of d{p^) is the same as the zero-bin contribution, M°, in calculating 
the integrated jet function in Eq. (41). The net collinear contribution from Fig. (b) is 
given as 


M, = M,- M° 


asCf 

277 


In 
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2 :euv 


1 1 / , 
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(62) 
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The naive collinear contribution from Fig. (c) is given as 


~ agCp,. 




271 

otxCj, 


— In 


47r 

QH'^ 


r(l — e) (p 2 )l+e Jp 2 /w'^ 

- 1 - 2 In 2 
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1 , fY 


{^(P^)(-— - In 
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2 In 


2QH 


QH^ 


+ 


r 1 QH^ - P 


M2 + ^p'^QH^+p'^ 

The zero-bin contribution is also suppressed, and is set to zero. 
The amplitude for the unintegrated jet function is given by 

/i) = 2(M, + M,) + Me + 4'^ + 


M 2 


}■ 


(63) 
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1 /3 


Le 


euv ^2 


+ In 




-- + 31n2- 
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(64) 


p^ \2 p2 _|_ Q2f2, 

If we put M = Qt, the last line vanishes. Dropping the UV divergence, the hnal unintegrated 
jet function at order as in the cone-type jet algorithm is given as 


UsCp (/3 , 




277 
r 2 

-p 2 


QH^ 
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In 
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QH^ 

1(3. 

V9 


+ b„= 


7 „, „ 

QH^ ^ 2 ^ 


3g2t2 


p 2 + QH‘^' 




(65) 


Note that the coefficient of 6{p‘^) is exactly the same as the integrated jet function, and 
the integration of the M^ distribution yields zero. Therefore we have successfully obtained 
the result that the integrated jet function is indeed obtained by integrating the unintegrated 
jet function. 


B. Unintegrated jet fnnction in the Sterman-Weinberg jet algorithm 


We can proceed to compute the unintegrated jet function in the SW jet algorithm with 
the same idea and method as in the cone-type jet algorithm. The virtual contributions from 
Fig. 1^ (a) and the self-energy correction are the same. The collinear contribution from Fig. 
(b) is given as 


M, = M,- M° 


ttsCp 
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( 66 ) 
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Here S'i(M^) is given by 


5i(M2) = Li2(w+) - Li2(M;-) - 2 J 1 - 


4M2 


w+ 


QH^ ^ w- 


, M 2 , M^/QH^ , 1,2 1,2 


where Li 2 (x) is the dilogarithmic function and 

1 


w± = 


{i±^ 


UP 


( 67 ) 


( 68 ) 


Later we will choose M = Qt/2, which corresponds to M'^= 1/4, w± = 1/2. And 

Si{QH‘^/A) = 0. 

Fig.g (c) is given as 

,2 Q 


M,= 


a.C 


s'^F 

2f 


{dip^ 


^ In 


2 eiR 2 QH^ 2 


+ ^2 (M 2 
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r 1 

-V 


4p2 


QH^ 
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where 


^ 2 (M 2 ) = Wl - 


4M2 1 , w+ 

- - In —, 


(69) 

(70) 


QH"^ 2 ’ 

with S 2 {Q‘^t^/^) = 0 . 

The hnal collinear matrix element for the unintegrated jet function in the SW algorithm 
is written as 

H + + . I ...2 

V 2 ^ QH^) ^ 2 


M~ = 

271 


a.Cpr., oM 1 1 /3 . /i^ X 3 


-{<^(P^) 
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'-^uv euv^2 


+ - hT 
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( 71 ) 


where S{M‘^) is the function which vanishes when M = Qt/2, and is given by 

S{M^) = 2Si{M^) + S2{M‘^) 

Q 

= 2 Li 2 (w+) - 2 Li 2 (w-) - 3^1 - AM'^/QH'^ + - In — 

V V 


W- 


+ In W- — In w+ — 2 In 


2 , M 2 M^/QH‘^ 

2„.. 0 1„ 


Q 2 t 2 w+-M^/QH‘^' 

The unintegrated jet function at order as in the SW jet algorithm is given by 


(72) 


asCp ( o. /3 , /i^ 
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1 + y^l — Ap'^/QH'^. 


)] (73) 

J\M=Qt/2\ ^ ' 


1 — \J^ — /QH^' 

Here the relation between the unintegrated and the integrated jet functions remains the 
same as in the cone-type algorithm. 
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VI. SOFT FUNCTION 


The soft function with a jet algorithm is defined as 

56 = ^ ^TV(o|yJn9.„ftK„'n|o>. 


(74) 


The Feynman diagrams for the soft function at order as are shown in Fig. Fig. (a) and 
(b) describe the virtual and real contributions respectively and the dashed line represents 
the cut. The soft function is twice the contributions of Fig. since the hermitian conjugate 
should be added. 

The virtual correction is given as 
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agCf 

2n 


poo 

(e'^®/i^)T(l + e) / dudviuv) 
Jo 


\-l-e 


a.gCp ^ 1 


271 '^euv eiR 


(75) 


The real gluon emission can be computed according to the jet algorithm in Eq. (29). The 


phase space available for the soft function is given in Fig. First, the contribution in which 
the soft gluon is in the jet in the n direction is given by 


Ml = 


agCp 

271 F(1 — e) 


dlJ 


-l-e 


/2/3Q 



dl+l-^-\ 


(76) 


Here the region /_ < 2j3Q is excluded since it will be included when we consider the case when 
the energy of the soft gluon is less than jdQ. Since the Z_ integral has the UV singularity, 
e > 0, the integral is rewritten as 




dl+l 
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(77) 


(78) 



FIG. 5. Feynman diagrams for soft functions at one loop (a) virtual corrections and (b) real gluon 
emission. The diagrams with the hermitian conjugate are omitted. 
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The amplitude for the case in which the soft gluon in the n jet is the same as Mi. 
The contribution from the soft gluon with energy less than f^Q is given as 

asCp , /■2/3Q-h 
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Summing over all the contributions, the real contribution is given by 


Sb = 2Mi + M 2 = 
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euv 2/3Q 2(3Qt 

There are two small regions missing in the above computation 2/3Q/(l +t^) < /+ < 2/3Q 

in Fig. but the contributions from those regions are negligible for small t ~ 0{Xy. 

The amplitude for the soft function is given by 
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M,oft = 2(^, + Sb) = 
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and the integrated soft function at order as is given as 
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VII. JET CROSS SECTIONS 


Let US summarize all the results for the factorized parts before we discuss the dijet cross 

Ttt^' 


sections. The hard function at order is given from Eq. (10) as 


«“'(Q".fi) = 


27r Q 2 Q 

The integrated jet functions at order a* are given by 
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The soft function at order as is given as 
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The anomalous dimensions of the hard, collinear and soft functions are given as 
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which yields 


'yu + 27j + 7s = 0. 


(87) 


It means that the dijet cross section is independent of the renormalization scale /r to next- 
to-leading log order. Explicitly at order the dijet cross section in the SW algorithm is 
given by |6] 

4w = <To(7f">(Q".A.) + 

- <7o-^^f-41n2/31ni - 31ni + ^ - 7)> (88) 

TT 8 / d 2 

whereas the dijet cross section in the cone-type algorithm differs only in the constant term 
which is —l/2-|-31n2 instead of 5/2 — 


VIII. DISCUSSION ON VARIOUS JET FUNCTIONS 


There have been many forms and definitions about the jet functions. Here we have 
confirmed that the integrated jet function is obtained by integrating the unintegrated jet 
function. On the other hand, the unmeasured jet function in Ref. m is not obtained by 
integrating the measured jet function. Especially, the UV divergence, or the corresponding 
anomalous dimensions are different in the measured and unmeasured jet functions. This may 
be confusing at first sight, but the definitions and the power counting are different, therefore 
there is no inconsistency. For example, when the total dijet cross section is computed, the 
calculations using the unintegrated, integrated jet functions, and the measured, unmeasured 
jet functions (with the corresponding soft functions) yield the same result. Here we discuss 
the origins of the difference and physics in detail. 

First, the inclusive jet function is the jet function without any jet algorithm and it has 
been computed to one loop [22H21], and to two loops |2S]- At one loop, the divergent part 
of the inclusive jet function is given as 
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^uv ‘^uv v 
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where only the divergent terms proportional to are shown. In this form, it is easy 

to compare the difference of the anomalous dimensions when integrated over the relevant 
variables. 

The corresponding divergent terms in the unintegrated jet function and the measured 
and unmeasured jet functions with the cone-type algorithm are given as 
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where tq = for a = 0 in Ref. [TT]. In comparing the differences, it is convenient to 

note the coefficients of the double poles of the various jet functions. 


The integrated jet function and the unmeasured jet function are the same including the 
hnite terms. Actually these dehnitions are the same. The difference lies in the unintegrated 
and the measured jet functions. In computing the unintegrated jet function, the jet algorithm 
requires that two particles be in the cone size of order f ~ A. In contrast, the measured 
jet function requires that the particles contributing to r should be conhned in a cone with 
the size of order \/t^ ~ A, while the cone size is approximately i? ~ A°. It means that in 
the measured jet function, there is effectively no jet algorithm involved, and the jet size is 
controlled only by the physical observable r^. The difference of the jet conhgurations in the 
unintegrated jet function and in the measured jet function is sketched in Fig. Therefore 
the idea of the inclusive jet function (with no jet algorithm) is employed in the measured jet 
function as far as the jet algorithm is concerned. When compared, it is clear from Eqs. (89) 
and (90) that the divergent parts of the inclusive jet function and the measured jet function 
are the same. This difference causes the different behavior of the divergences between the 
unintegrated and the measured jet functions. 


Then there also arises the question of the behavior of the dijet cross section,which is a 
convolution of the hard coefficient (the same), the jet functions, and the soft function. When 
the integrated or the unintegrated jet function is used, since their anomalous dimensions 
are the same, the dijet cross section is the same, and independent of the renormalization 
scale /i. However, it seems that the cross sections with the unmeasured jet function and the 
measured jet function integrated over r may have different anomalous dimensions, which is 
troublesome. The answer lies in the fact that the measured soft function is also affected by 
the jet size of a/t ~ A with t ~ A*^. If the change of the measured soft function is included, 
we have verihed that the same dijet cross section is obtained. Therefore the dijet cross 



FIG. 6. Comparison of the jet sizes according to the power counting of the parameters, (a) In 
the unintegrated jet function, the cone size (solid line) is of order t ~ A. (b) In the measured jet 
function, the jet cone size is of order i? ~ A^, while the physical cone size (dashed line) is controlled 
by Vt A, which is of the same order as t in (a). 
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section remains the same whichever jet functions with their corresponding soft functions are 
used. 


IX. CONCLUSION 

The implementation of the jet algorithms in SCET successfully factorizes the dijet cross 
section in e’''e“ annihilation. We have confirmed that the divergences in each factorized 
part are truly of the UV origin in the cone-type and the SW jet algorithms. Thus the 
jet function and the soft function are infrared hnite in each jet algorithm. This is verihed 
in pure dimensional regularization and in the regularization with rapidity regulator in the 
collinear Wilson line. The question still remains whether it is also true in various other jet 
algorithms such as the exclusive [H] and Georgi [22] algorithms. It will be presented soon 
elsewhere m 

In proving that the jet and the soft functions are IR hnite in pure dimensional regulariza¬ 
tion with the spacetime dimension D = 4 — 2e, the IR and UV divergences are entangled due 
to the phase space constraint. The phase space constraints for the collinear, zero-bin and 
soft contributions complicate the identihcation of the sources of the divergence. We have 
performed a careful analysis to identify the source of the divergences and to disentangle the 
IR and UV divergences. And in the cone-type and in the SW jet algorithms, the jet function 
and the soft function are proven to be IR hnite to next-to-leading order, which enables the 
factorization of the dijet cross section. 

The integrated jet function is indeed obtained by integrating the unintegrated jet function. 
In this case, as was discussed in comparing with the relation between the measured and 
unmeasured jet functions, the power counting of the parameters involved is important. The 
power counting for the (un)integrated jet function is such that f ~ A, which is related to the 
cone size, and {3 ~ A^, which is related to the jet veto. 

In calculating the unintegrated jet function, it is important to keep the power counting 
correctly. Especially in computing the zero-bin contribution M°, the distribution function 
should be concentrated near = 0 by power counting. This distinction makes the relation 
between the unintegrated and the integrated jet functions correct. If we do not care about 
the power counting and compute the zero-bin contribution with the distribution away from 

= 0, can have any value and the result becomes exactly the inclusive jet function with 
t = I, instead of the unintegrated jet function. This is obvious because the inclusive jet 
function includes all the regions of p^. However, if we employ the inclusive jet function, 
or if we neglect the careful power counting, the dijet cross section would depend on the 
renormalization scale p, which is not sensible. Therefore careful power counting should be 
performed in computing jet functions. 

It will be interesting to probe the divergence structure in other jet algorithms not only in 
e’''e“ annihilation but also in hadronic collisions. The kinematics in both cases is different, 
and the color structure is more complicated since the initial states also consist of hadrons. 
The factorization property with various jet algorithms in different scattering processes will 
be pursued based on the analysis in this paper. 
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Appendix: Integrated jet functions in the SW jet algorithm with rapidity regulator 


Jet functions can be computed with a rapidity regulator for the IR divergence, while the 
UV divergence is handled by the dimensional regularization. For the integrated jet function, 
we use Eq. ( [3^ in the collinear Wilson line to regulate the IR divergence. The virtual 
corrections in Fig. (a) depend on the regulator A, but the net collinear contribution after 
the zero-bin subtraction yields the same result as in Eq. (37). 

In the SW jet algorithm, the naive collinear contribution in Fig. (b) is given as 
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For the zero-bin contribution, since the IR divergence is handled by the IR regulator, the e 
should be of the UV origin, e > 0. Therefore it yields 
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And the net collinear contribution is given by 
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Likewise, M^ is given by 
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Finally the total collinear amplitude is given as 
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which is the same as Eq. (50). 
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